Self— Consistent Conserving Theory for Quantum 
Impurity Systems: 
Renormalization Group Analysis* 

Stefan Kirchner and Johann Kroha 



Institut fur Theorie der Kondensierten Materie, Universitdt Karlsruhe 
Postfach 6980, 76128 Karlsruhe, Germany 



We review the diagrammatic, conserving theory for quantum impurities with 
strong on-site repulsion. The method is based on auxiliary particle technique, 
where Wick 's theorem is valid, which opens up the possibility for generaliza- 
tions to more complicated situations. An analysis in terms of the perturbative 
renormalization group (RG) shows that on the level of the Conserving T- 
matrix Approximation the theory correctly describes the RG flow, including 
the non-scaling of potential scattering terms and the correct Kondo temper- 
ature. 

PACS numbers: 11.21. +a, 71.10.Fd, 75.20.Hr, 11. 10. Hi. 



1. INTRODUCTION 

The physics of correlated electron systems on a lattice, such as heavy 
fermion compounds or narrow band metals is determined by a strong on- 
site repulsion between the electrons. It can lead to the formation of local 
magnetic moments and, subsequently, either to the Kondo effectu or to mag- 
netic ordering, to a Mott-Hubbard metal-insulator transition, and even to 
superconductivity. The theoretical description of such systems has been 
greatly advanced by the dynamical mean field theory (DMFT)J3tJ where 
the lattice system is mapped onto an effective single-impurity Anderson 
model, embedded self-consistently in a correlated electron bath. Thus, the 
complexity of the lattice system is reduced to that of a quantum impu- 
rity problem. Quantum impurity systems are also of central interest in 
their own right: The Anderson impurity model is the standard model, e.g., 
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for the description of quantum dots in the Coulomb blockade and Kondo 
regimesElli, of scanning tunneling spectroscopy (STS) on magnetic ions on 
a metal surface^™, and for the investigation of single-impurity physics in 
heavy fermion systemsclJO. Therefore, it is highly desirable to construct 
theoretical methods flexible enough to describe quantum impurities with a 
complex structure of the fermionic bath (in DMFT), with multiple local or- 
bitals (in rare earth systems, quantum dots etc.), and which can readily be 
generalized to non-equilibrium transport. Existing exact solution methods 
like Bethe ansatz, conformal field theory (CFT), and numerical renormal- 
ization group (NRG) are lacking this flexibility to a large extent. 

A general method, based on the auxilary particle technique^! (where 
Wick's theorem is valid), which is capable of describing quantum impurity 
systems over the complete temperature range from the weak to the strong 
coupling regime and which, as a diagrammatic method, offers the above- 
mentioned flexibility, has been envisaged by Peter Wolfle starting in the late 

^al break- 
possible. 



1980s, even before the advent of DMFT and before the techno 
through that made the fabrication of Kondo quantum dotsH3' 
Such a method, now termed Consuming T-Matrix Approximation (CTMA), 
has subsequently been developeduO'EJ'Ej, and has been shown to describe, 
as an example of physical quantities, the spin susceptibility of the single- and 
the two-channel Anderson impurity model correctly from the high temper- 
ature regime down to thejowest temperatures T considered of about 1/100 
of the Kondo scale Tk e9>EH. 

The choice of diagramsmmprising the CTMA has been justified in 
terms of principal diagramsEJ'", where in each order of the impurity hy- 
bridization the dominant term in the spin as well as in the charge fluctua- 
tion channel are taken into account. In the present paper we present another 
check for the validity of the CTMA: The renormalization group flow of the 
coupling constants appearing in the spin and charge vertices defined by the 
CTMA should be the same as the exact flow of the respective coupling con- 
stants of the Anderson impurity model at least in the weak coupling regime. 
By self-consistency, the CTMA may then be expected to reach the correct 
strong coupling behavior. We note that this scheme might be extended to 
be used as a guiding principle for the construction of approximative strong 
coupling methods for other models like the Hubbard and t-J-modelscJ. 

The paper is organized as follows. After a brief overview of the auxil- 
iary particle technique and its exact properties in Section 2, the conserving 
approximations for the Anderson impurity model, the Non-Crossing Ap- 
proximation (NCA) and the Conserving T-Matrix Approximation (CTMA), 
will be defined in Section 3. Section 4 contains the perturbative renormaliza- 
tion group analysis of both the NCA and the CTMA in the sense discussed 
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above. Some conlusions are drawn in Section 5. 



2. EXACT PROPERTIES OF THE AUXILIARY PARTICLE 

TECHNIQUE 

The large on-site repulsion U between electrons in the local d-orbital of 
an Anderson impurity effectively restricts the dynamics to the sector of Fock 
space with no double occupancy. It can be implemented using the auxiliary 
or slave boson methodEl, where in the limit U — » oo the creation operator 
for an electron with spin a in the cf-level is written in terms of the auxiliary 
fermion operators f a and boson operators b as S a = f^b. This representation 
is exact, if the constraint on the total auxiliary particle number operator, 
Q = J2a fafa + = 1, is fulfilled, fl and ¥ may be envisaged as creating 
the three allowed states of the impurity: singly occupied with spin a or 
empty. The Hamiltonian of a single Anderson impurity with local level Ed, 
embedded in a sea of conduction electrons (creation operators ct , dispersion 

kcr 

er, half band width D, and density of states at the Fermi level Nq) via a 
hybridization matrix element V then reads (U — > oo), 

k,cr a k,cr 

In the Kondo regime, A^o Jo = NoV 2 /\Ed\ <C 1, the model has a low temper- 
ature scale, the Kondo temperature Tk — De~ l l( 2NoJo \ at which the system 
crosses over from singular spin scattering to a spin-screened, strong coupling 
Fermi liquid (FL) ground state.El 



2.1. Exact Projection onto the Physical Fock Space 

The auxiliary particle Hamiltonian (|l]) is invariant under simultaneous, 
local U{1) gauge transformations, f a — * fa^^, b —> be % ^ T \ with </>(r) 
an arbitrary, time dependent phase. While the gauge symmetry guarantees 
the conservation of the local, integer charge Q, it does not single out any 
particular Q sector, like Q = 1. In order to projectonto the Q = 1 sector 
of Fock space, one may use the following procedureEfrEil: Consider first the 
grand-canonical ensemble with respect to Q and the associated chemical 
potential —A. The expectation value in the Q = 1 subspace of any physical 
operator A acting on the impurity states is then obtained exactly as 
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where the index G denotes the grand canonical ensemble, Q denotes the fu- 
gacity £ = e~^ A , and — A is by construction the chemical potential associated 
with the local charge Q. In the second equality of Eq. @ we have used the 
fact that any physical operator A acting on the impurity states is composed 
of the impurity electron operators d a , S a , and thus annihilates the states 
in the Q = sector, A\Q = 0} = 0. It is obvious that the grand-canonical 
expectation values involved in Eq. (||) may be factorized into auxiliary par- 
ticle propagators using Wick's theorem, thus allowing for the application of 



standard diagrammatic techniques. For a detailed review see Ref. 21,p2|. 

It is important to note that, in general, A plays the role of a time 
dependent gauge field. In Eq. (||) a time independent gauge for A has been 
chosen. In this way, the projection is only performed at one instant of time, 
explicitly exploiting the conservation of the local charge Q. This means 
that in the subsequent development of the theory, the Q conservation must 
be implemented exactly. It is achieved in a systematic way by means of 
conserving approximations^], i.e. by deriving all self-energies and vertices 
by functional derivation from one common Luttinger-Ward functional <5 of 
the fully renormalized Green's functions, 

S 6)/)C = 5${G b , G f , G c }/5G bJjC . (3) 

This amounts to calculating all quantities of the theory in a self-consistent 
way, but has the great advantage that gauge field fluctuations need not be 
considered. 



2.2. Infrared Threshold Behavior of Auxilary Propagators 

As seen from Eq. (|2|), the limit A — > oo effecting the projection onto 
the physical subspace implies that the traces involved in the time ordered 
pseudofermion and slave boson Green's functions Gf,Gb are extended purely 
over the Q = sector of Fock space, and thus the backward-in-time contribu- 
tion to the auxiliary particle propagators vanishes. Consequently, the auxil- 
iary particle propagators are formally identical to the core hole propagators 
appearing in the well-known X-ray problem@, and the long-time behavior of 
Gf (Gb) is determined by the orthogonality catastrophe^ of the overlap of 
the Fermi sea without impurity (Q = 0) and the fully interacting conduction 
electron sea in the presence of a pseudofermion (slave boson) (Q = 1). It 
may be shown that the auxiliary particle spectral functions have threshold 
behavior with vanishing spectral weight at T = for energies u below a 
threshold E a , and power law behavior above E Q , At b {oj) oc Q(lo — E )cu~ a f' b . 

In the present paper we will consider only the single-channel Anderson 
model with its spin screened FL ground state, leaving the impurity as a 
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pure potential scatterer. In this case, the exact threshold exponents may be 
deduced EUlZEl from an analysis in terms of scattering phase shifts0, using 
the Friedel sum rule. One obtains for spin degeneracy N > 1 and conduction 
channel degeneracy M = 1, 

a f = ^p^, « 6 = l-f (N>1,M = 1), (4) 

These results have been confirmed by numerical renormalization group (NRG) 
calculations^ and by use of the Bethe ansatz solution in connection with 
boundary CFTeHI. We note in passing that, on the contrary, in the non- 
Fermi liquid (NFL) case of the multi-channel Kondo model (M > N > 2) 
the threshold exponents have been deduced by a CFT solutionis as 

a { = ^ M Ar , a b = — - A7 (M > N > 2) . (5) 

} M + N ' M + N \ - - j \ i 

Since the dependence of a/, a& on the impurity occupation number rid shown 
above originates from pure potential scattering, it is characteristic for the FL 
case. The auxiliary particle threshold exponents are, therefore, indicators 
for FL or NFL behavior in quantum impurity models of the Anderson type. 



3. CONSERVING APPROXIMATIONS 

The task in constructing a self-consistent, conserving theory for a given 
system consists in finding the proper Luttinger-Ward functional that yields 
a correct description for the respective strongly correlated model. In this 
section we describe two approximation schemes for the Anderson impurity 
model, the NCA and the CTMA. 

In Matsubara representation the auxiliary particle Green's functions 
Gf : b read in terms of the self-energies S/^, 

Gf, b (iuj n ) = {[Gy ^iu^r 1 - Ef }b (iu n )\ , (6) 

where 

G%(iiu n ) = (iu n -E d - A)" 1 , G° b (iu n ) = {ioj n - A)" 1 (7) 

are the respective free propagators. Since, as a consequence of the projection 
procedure A — > oo, the energy eigenvalues of H + XQ scale to infinity cx XQ, 
it is useful to take a gauge A = Ao + A' with A' — > oo, to shift the zero of 
the auxiliary particle frequency scale by A', and to determine Ao such that 
for T = the threshold of the auxiliary particle spectral functions lies at 
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Fig. 1. Diagrammatic representation of the generating functional of 
the NCA. Also shown are the pseudoparticle self-energies and the local d- 
electron Green's function derived from <£, Eqs. (11)— (13). Throughout this 
article, dashed, wavy and solid lines represent fermion, boson, and conduc- 
tion electron lines, respectively. 

the frequency uj = 0.0 In particular, for vanishing hybridization the Green's 
functions read in this gauge, 

G%(ioj n ) = (iun)- 1 , Gg(iw„) = (iw„ + EdY 1 . (8) 

3.1. Non Crossing Approximation (NCA) 

The NCA as the simplest conserving approximation for the Anderson 
impurity model is motivated as a self-consistent expansion in terms of the 
hybridization parameter V which is assumed to be small. The NCA gen- 
erating functional <J> is shown in Fig. |l| Also shown are the self-energies 

and the physical Green's function for an electron in the local d-level, 
Gd, generated from this by functional differentiation. They obey after an- 
alytic continuation to real frequencies (ico — ► to — iO) and projection onto 
the physical subspace the following equations of self-consistent second order 
perturbation theory 

T / ^f(e)A° ca (-e)G b (u + e-iO) (9) 
r£ / -f(e)J%,(e)G f< ,(u, + e-iD) (10) 

J de e~ PE [G ftr (w + e - iO)A b (e) 

-A fa (e)G b (e-iv + iO)} , (11) 

where T = vrA^o^ 2 , and A® ca = -ImG^/No is the (unrenormalized) con- 
duction electron density of states per spin, normalized to the density of 
states at the Fermi level Nq, Af^{oS) = hnGf t b(u) — iO) denote the imagi- 
nary parts of the advanced propagators, and /(e) = l/(exp(/3e) + 1) is the 
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Fig. 2. Diagrammatic representation of the Bethe-Salpeter equation for 
(1) the conduction electron-pseudofermion p-h T- matrix T^f) A , (2) the con- 
duction electron-pseudofermion p-p T-matrix JW', Eq. (|l2|), and (3) the 
conduction electron-slave boson T-matrix T"( cfc ) , Eq. (|i~3|) . The external lines 
are drawn for clarity and do not belong to the T-matrices. 



Fermi distribution function. Together with the expressions (||), (||) for the 
Green's functions, Eqs. (|9|)~(|Tl|) form a set of self-consistent equations for 
S&,/ c > comprised of all diagrams without any crossing propagator linesllil. 
For an efficient procedure for the numerical evaluation of the self-consistent 
equations at low T see Ref. The solutions of the NCA equations have 
threshold power law behavior, with the exponents given by Eq. (H|),0 which 
are in disagreement with the correct exponents for the FL case, Eq. (||). As 
a consequence, the NCA does not correctly describe the FL regime of the 
Anderson impurity model, producing a spurious low energy singularity in 
the local (i-electron spectrum, although it gives qualitatively correct results 
at high T and down to temperatures of about the Kondo scale Tr-. 



3.2. Conserving T-Matrix Approximation (CTMA) 

In order to construct an aproximation which eliminates the shortcom- 
ings of the NCA mentioned above, the guiding principle has been to include 
those contributions to the vertex functions which capture the FL physics of 
the single-channel Anderson model below i.e. the formation of a collec- 
tive singlet bound state of the local and of the conduction electron spins J 1 ^r 1 ! 
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spin 
fluctuations 
charge 



Fig. 3. Diagrammatic representation of the Luttinger-Ward functional 
generating the CTMA. The terms with the conduction electron lines running 
clockwise (labelled "spin fluctuations") generate while the terms with 

the conduction electron lines running counter-clockwise (labelled "charge 
fluctuations") generate T^ cb \ The two-loop diagram is excluded, because it 
is not a skeleton. 



The onset of a bound state is usually evidenced by a pole in the correspond- 
ing two-particle correlation function or T-matrix. Hence, one expects a pole 
in the singlet channel of the conduction electron-pseudofermion T-matrix 
Taking a single slave boson propagator as the particle-particle (p-p) 
and particle-hole (p— h) irreducible c — f vertex (leading order in V con- 
tribution), the equations for the c — / p-h and p-p T-matrices are shown 
diagrammatically in Fig. [2] (1), (2). These contributions include, at any given 
order in the hybridization V, the maximum number of spin flips possible. 
Considering that in the mixed valence regime the Anderson model is domi- 
nated equally by spin and by charge fluctuations, one should also take the 
conduction electron-slave boson T-matrix T^ cb ' into account, which simi- 
larly includes the maximum number of charge fluctuation processes (see Fig. 
^ (3)). The (linear) Bethe-Salpeter equations for T^ c f \ read, 

T^fl, T ,(iu n ,iLo' n , iQ n ) = + V 2 G b {iuj n + iJ n - i£l n )5 ar i5 TCT i 

- V 2 TY,G b {iLo n + iJ^-i^l n )x (12) 

< 

° {iuj n , iJ n , in n ) = + v 2 G fa {+iuo n + %J n - in n ) 

- V^Y^Gf^n + iJ^-iVLn) x (13) 

g 6 («4) g° ct (-^" - »«n) r (cb) °(iu';M n ,in n ). 

Inserting the NCA solutions for G j a , G b into Eq. (|l2|), one finds after an- 
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Fig. 4. Diagrammatic representation of the CTMA pseudofermion and slave 
boson self-energies £/, 



alytic continuation, as expected, a pole in the singlet channel of T^ c ^ at 
a center-of-mass frequency £Iq which scales with the Kondo temperature, 
Q Ri —Tk- In order to incorporate this physics into the self-consistent 
scheme, we must construct a Luttinger-Ward functinal $ which includes 
the total vertices T^^ A , T^ c ^\ on the level of the self-energies. This 
functional is shown in Fig. |3| and consists of the sum of all slave particle rings 
where the bare hybridization vertices are connected by conduction electron 
lines in such a way that at most two other hybridization vertices are spanned. 
The set of self-energy terms deduced from by functional differentiation 
is shown in Fig. ^|. In the self-consistent scheme all propagators appearing 
in the self-energies are understood to be the fully renormalized ones. In 
order to avoid double counting of terms, the one-rung term in the second 
diagram of £ j and the one- and two-rung terms in the third diagram of £ f , 
Fig. must be excluded. An analogous exclusion of terms must be done in 

p— i ( c f) 

Sb. For a detailed description see Ref. |21]. Note that T^ CA just corresponds 
to a renormalization of the boson propagator by the NCA self-energy and, 
hence, is already included on the NCA level of self-consistent approxima- 
tion. The self-consistent set of non-linear integral equations resulting from 
the definition of the self-energies Fig. |] and the full Green's functions Eq. (JgJ) 
is called CTMA. Solving the CTMA equations selfconsistently by iteration, 
the aforementioned pole is shifted to the threshold frequency $7 = where 
it merges with the continuous spectral density and thus renormalizes the 
threshold exponents. As signatures of the correct FL behavior below Tk it 
was found from the numerical solutions of the CTMA equations first that the 
correct FL exponents of the auxiliary particle propagators, Eq. (||), are re- 
produced within the error bars of the numerical solution, and second that the 
static impurity spin susceptibility x shows Pauli behavior (x(T) ~ const.) 
for T < Tk down to the lowest T considered, T « 10~ 2 • Tk (see Fig. ||). 
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Fig. 5. Left panel: The fermion and boson threshold exponents a/, ol\, are 
shown for N = 2, M = 1 in dependence of the average impurity occupation 
rid- Solid lines: exact values, Eq. (||); Symbols with error bars: CTMA; 
dashed lines: NCA. Right panel: Static susceptibility of the single-channel 
Anderson impurity model in the Kondo regime (E^ = —0.8D, T = 0.1D, 
Lande factor g = 2). The CTMA result is calculated using the The CTMA 
and NCA results are compared to the Bethe ansatz result for the Kondo 
modelEa. The CTMA susceptibility obeys scaling behavior in accordance 
with the exact results (not shown). 




4. PERTURB ATIVE RG ANALYSIS 

Having introduced the CTMA as a method to describe quantum im- 
purities in the strong coupling region, it is also important to know how it 
performs in the weak to intermediate coupling region, i.e. how within the 
CTMA the coupling constants for spin and potential scattering are renor- 
malized in the sense of the perturbative renormalization group (RG). To 
obtain the correct weak to intermediate coupling renormalization within a 
self-consistent treatment turns out to be a non-trivial task. As seen below, 
not even within the NCA, which is a self-consistent expansion in the small 
parameter V (hybridization), the correct weak coupling renormalization is 
produced for the Anderson impurity model. 

In this section the coupling constant renormalization for spin and po- 
tential scattering will be analyzed within the perturbative RG0cl. In doing 
so we will assume to be in the Kondo regime of the Anderson model, where 
E d < and F/\E d \ <C 1. 

In the (anisotropic) Kondo model the local spin-conduction electron 
interaction vertex is defined as 

J±(S + a- + S-a + ) + J\\S z a z . (14) 

For isotropic coupling, J± = Jn = J, this reduces to the regular Heisenberg 
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(2) 



Fig. 6. (1) The fermion selfenergy expressed in terms of the exact 3-point 
vertex V (black triangle) is shown. (2) Representation of V in terms of the 



4-point vertices 7 C ^, 7 cfe 
bare 3-point vertex V. 



The first diagram on the right-hand side is the 



coupling J S ■ a. In the above expressions S = (S x , S y , S z ) denotes the 
impurity spin 1/2 vector operator and a = (cr x ,ay,a z ) the vector of Pauli 
matrices. S and are the spin raising (+) and lowering (-) operators 
for the local and for the conduction electron spin, respectively. In addition, 
there may be a potential scattering term which has the structure 

W(ll) (15) 

with 1 the 2x2 unit matrix in impurity and in conduction electron spin 
space. 

Since in the Anderson model the hybridization V (Eq. ([l])) appears 
as the bare coupling constant, but not the spin and potential scattering 
couplings Jj_, Jm, and W, the NCA or CTMA self-energies must first be 
represented in terms of the 4-point spin and potential vertices. This is done 
as follows for the pseudofermion self-energy Sj. We note that it is not 
necessary to consider separately the boson self-energy £{, here, since will 
appear naturally as a coupling constant renormalization in the perturbative 
RG. The (exact) £/ can be written in terms of the exact, total 3-point 
hybridization vertex V in slave boson representation as (see Fig. || (1)) 

£ /(r (w-iO) = V J ^f(e)G° b (e + u-iO) x (16) 
[V(u - iO, e + iO) G° ca (-e - iO) - V{uo - iO, e - iO) G° ca (-e + iO)] , 

where G® is the bare boson Green's function, Eq. (@). The total 3-point 
vertex V can in turn be expressed in terms of the exact, total 4-point c-f 
and c-b vertices j c f , j cb , as shown in Fig. || (2). These 4-point vertices 
are comprised of all vertex corrections connecting the f and the c line or 
the b and the c line of the 3-point vertex V, respectively. Note that vertex 
corrections connecting the f and b lines do not appear as these are of higher 
order in the fugacity e _/3A and thus vanish under projection onto the physical 
subspace. A given approximation for the selfenergy S / (like NCA or CTMA) 
corresponds to a respective approximation for 7 C ^, j cb . 
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(1) (2) (3) 

Fig. 7. perturbative RG renormalizations of the (irreducible) c-f vertex. 
Diagram (1) is the NCA result, the sum of diagrams (1), (2) the result of 
CTMA. (3) ist the contribution from the c-b vertex, which is not scaling. 
The black dots on the conduction electron lines indicate that the frequency 
integrals in those lines are restricted to the region [-D, —D + dD]. 



The perturbative RG analysis can now be applied to the Bethe-Salpeter 
equations for the total vertices (T-matrices) 7 C ^, j cb , thus determining the 
RG flow of the spin and potential scattering coupling constants within the 
given approximation. Note that for the perturbative RG treatment all propa- 
gator lines appearing are the bare ones, while in the self-consistent treatment 
described above all propagators are fully renormalized. 



4.1. Perturbative RG for the NCA 

(cf) f 

Inserting the "RPA-like" T-matrix T^ CA for 7 C/ and the resulting V 
into Eq. ( |l6|) for the f self-energy, it is seen that the NCA is just defined by 

7 C/ = T^ c \ (17) 
l cb = 0, (18) 

(cf) 

and 2jvca P rov ides just the renormalization of the boson propagator by the 
NCA self-energy, as mentiond above. Integrating out the high energy region 

(cf) 

from the Bethe-Salpeter equation JnL.T^ A and applying the usual pertur- 
bative poor man's scaling anarysiseilc3, we obtain for the renormalization of 
the irreducible conduction electron-pseudofermion vertex under a change of 
the high energy cutoff from D to D — dD (compare Fig. 0(1)) 

dA (cf) = _ No ^R £ [A^J^^IA^]^^ (19) 

s"o" 

where the spin indices s, s', s", a, a f , a", are as defined in Fig. (1) and 
the (unrenormalized) c-f vertex is given by the inhomogeneous part of the 
Bethe-Salpeter equation (|l^) [see also Fig. | (2)], 

A [ fJ,sa = [JAS + °- + S~a+) + J n S z a z + Wll] sla ,, S(T . (20) 
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In Eq. ( j20|) s, s' denote the conduction electron spin indices, a, a' the impu- 
rity spin indices, where an index with (without) prime represents an ingoing 
(outgoing) particle. For the Anderson model we have J± = Ju = 2W = 
V 2 /\Ed\, the l-Erfl -1 arising from the bare boson propagator, Eq. ||, taken 
at frequencies \u\ <C \E<i\. However, for generality we keep the anisotropic 
couplings Jj_, Ju, W in the above equations. dA^ c ^ may be simplified as 

dAW = -N dlnD \ (J,, J ± + 2WJ ± )[S + a- + S-<r + ] 

+ {Jl + 2WJ\\)[S z a z ) (21) 



+ C^Jl + \j\\ + W 2 )[11] 



Hence, within NCA we have the RG equations 
dJ ± 



-N {J\\J± + 2WJ ± ) (22) 
dhxD -No(Jl+2WJ ll ) (23) 

- -m\* +l J +w *>- (24) 



dlnD 
dJu 



For the initial conditions of the Anderson impurity case, J± t o = Jiiq = 
2 Wo = F 2 /|Srf| = Jo, these are easily integrated to give 

J{D) = D • ( 25 ) 

1 + 2A^ J ln^ 

i.e. the spin coupling constant diverges at the Kondo temperature Tx = 
-Doe _1// ^ 2Ar ° Jo - ) . However, the potential scattering coupling is also renormal- 
ized under the RG flow. Thus, within NCA there is a spurious divergence of 
the potential scattering term at Tk as well. The fact that within NCA po- 
tential scattering is incorrectly treated on the same footing as spin scattering 
may be traced back to be the origin why (1) within NCA the asymmetry of 
the Kondo resonance comes out too large and why (2) NCA gives a qualita- 
tively wrong descripting of the Kondo resonance in a magnetic field E3. 



4.2. Perturbative RG for the CTMA 



We now consider the coupling constant renormalization under the RG 
flow within CTMA. The CTMA f self-energy (Fig. |j) is generated by insert- 
ing, in addition to the NCA contribution, the ladder T-matrices, Eqs. (12), 
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fll3|), for "f( c f\ j( cb \ This is seen by comparing Fig. || with Fig. i.e. we 
have in addition to the vertices of the previous section, 

7 C/ = (26) 
7 cfe = . (27) 

The resulting additional c-f vertex renormalization under cutoff reduction is 
(see Fig. @ (2)) 

dA^ = +iV ^ £ [A^W, S , V ,[A^W„, SCT , (28) 

s"cr" 

which can again be simplified as 

dA< c/ > = +N dlnD | (-J||J± + 2WJ ± )[S + a- + S _ <t + ] 

+ (-Ji + 2W,|)[5 Z( 7 Z ] (29) 

+ (^I + ^||+^ 2 )[H] 

Finally, we also need to investigate the renormalization due to 7^ (Fig. 
^ (3)). This term is obviously diagonal in spin space (potential scattering). 
Its amplitude is 

r-D+dD 

dA^ = -N V* / de x 

J-D 

1 (30) 



G%{u + e)G%{u' + £ )G° b (e-n) 
V A dD 



Q=0,<j,tj' ,e<.D<\E d \ 



or, with W = V 2 /{2\E d \), 

dW = -4N \E d \W 2 ^ . (31) 

This correction is not logarithmic (non-scaling), and the integration of Eq. 
(|3l| ) yields a nondivergent result. Hence, it need not be considered in the 
RG flow. 

Adding up the contributions Eqs. (|2l|) and (p9|) it is seen that the po- 
tential terms always cancel each other. Thus, there is no renormalization of 
the potential term. The resulting RG equations of the CTMA are 
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-2N Jl 



= . 



(33) 



(34) 



These are identical to the perturbative RG equations of the original Kondo 
model. This proves that the CTMA incorporates the complete Kondo physics 
also in the weak and intermediate coupling regime, where the perturbative 
RG is valid. 



We have given an overview of a diagrammatic, conserving approxima- 
tion the CTMA, designed to describe quantum impurity systems of the An- 
derson type in the strong coupling regime. By means of a perturbative 
renormalization group analysis it was shown that, in contrast to an earlier 
approximation (NCA) it recovers the correct RG flow in the weak to inter- 
mediate coupling region. In particular, it describes correctly, that potential 
scattering terms do not scale in the Kondo and Anderson problems, while 
the correct Kondo scale for the divergence of the spin dependent part is re- 
covered. Together with the fact that in numerical solutions of the CTMA 
the auxiliary particle threshold exponents and the spin susceptibility show 
the correct Fermi liquid behavior of the strong coupling regime, this pro- 
vides evidence that the CTMA correctly describes the Anderson model both 
in the strong and in the weak coupling regime on the same footing. It has 
also the potential to be generalized to describe more complicated situations 
and lattice problems by use of the DMFT. 
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